Interaction strengths in cuprates from the inelastic neutron-scattering measurements 
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The t — U — V — J model is used to describe the positions of the experimental peaks associated 
with commensurate and incommensurate structure of the magnetic susceptibility probed by neutron 
scattering in cuprate compounds. Assuming that the tight-binding form of the mean-field electron 
energy and the maximum gap can be obtained by fitting the angle-resolved photoemission spec- 
troscopy (ARPES) data, we have determined the strengths of the on-site repulsive interaction U, 
the spin-independent attractive interaction V and the spin-dependent antiferromagnetic interaction 
J from the positions of the commensurate and incommensurate peaks. 
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Introduction. It is evident from many 

experiments^ that the spin excitation spectrum 
in the superconducting state of cuprate compounds 
for low temperature and low energy consists of in- 
commensurate (IC) magnetic peaks at the quartet of 
wavevectors Q s — (tt(1± 5), 7r(l± 6)) where 8 represents 
the degree of incommensurability (the lattice parameter 
a = 1 and h — 1 are set to unity, Q 5 corresponds 
to (0.5(1 ± <5), 0.5(1 ± 5)) in reciprocal lattice units, 
r.l.u., and the total number of sites is N). S decreases 
with increasing the energy transfer uj and vanishes 
at energy uj res . This commensurate peak is called 
a magnetic resonance peak and it is centered at the 
antiferromagnetic wave vector = (7r,7r). It was 

suggested^— that the existence of the resonance peak 
in Bi2212 samples is related to two strong interactions: 
an on-site repulsion U which drives the system close 
to an antiferromagnetic instability, and a short-range 
antiferromagnetic interaction J related to the fact that 
cuprates are oxides of copper doped with various other 
atoms (doped antiferromagnetic Mott insulators). Since 
d-wave superconductivity on the cuprate square lattice 
does not rule out the existence of a spin-independent 
attractive interaction, we naturally arrive at the idea 
that the t — U — V — J model could be a possible 
theoretical scenario which fits together three major parts 
of the high-T c superconductivity puzzle of the cuprate 
compounds: (i) it describes the opening of a d-wave 
pairing gap, (ii) it is consistent with the fact that the ba- 
sic pairing mechanism arises from the antiferromagnetic 
exchange correlations, and (iii) it takes into account the 
charge fluctuations associated with double occupancy of 
a site which play an essential role in doped systems. 

In the one-layer approximation, the Hamiltonian of the 
two-dimensional t — J — U — V model is: 



ti.f 1t (l) 



where (i is the chemical potential. The Fermi operator 
Yia (V^'.ct) creates (destroys) a fermion on the lattice site 
i with spin projection a =t, J, along a specified direction, 
and n,i. a = ipj o-ipi.a is the density operator on site i with 



nearest-neighbor sites. The spin operator is defined by 
fi = ip\ a ~cf aai^i^i /2, where 77" is the vector formed by 
the Pauli spin matrices {a x , a y , <j z ). The terms in ([1]) 
represent the hopping of electrons between sites of the 
lattice, their on-site repulsive interaction (U > 0), the 
attractive interaction between electrons on different sites 
of the lattice (V > 0) and the spin-dependent Heisenberg 
near-neighbor interactions, respectively. 

The gap equation in the case of d-wave pairing 
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where V^p 



(2) 



and 



2V + 3 J/2, E(k) = yej 

Ak = A(cosfc x — cosfc l/ )/2 provides a relationship 
between the strengths of the V and J interactions. The 
mean-field electron energy has a tight-binding form 
f?k = h (cos k x + cos k y ) /2 + t 2 cos k x cos k y +t 3 (cos 2k x + 
cos2k y )/2 + t\ (cos 2k x cosk y + 0082^ cos k x )/2 + 
t$ cos 2k x cos 2k y — /i obtained by fitting the ARPES 
data with a chemical potential fi and hopping ampli- 
tudes ti for first to fifth nearest neighbors on a square 
lattice. Since the spin susceptibility and the two-particle 
Green's function share common poles, we can adjust 
the parameters U, V and J in such a way that the 
resonance peak solution of the Bethe-Salpeter 

(BS) equations for the spin collective mode of the 
Hamiltonian (JTJ at a wave vector Q^. The observed 
IC peaks are also poles of the spin susceptibility (or 
solutions of the BS equations), and therefore, we can use 
any of the IC resonances to obtain the exact strengths 
of the interactions. The calculated strengths could be 
tested using the positions of the other IC peaks observed 
on the same sample but at different transfer energy. 

Bethe-Salpeter equations for the collec- 
tive modes. Thought the method for solv- 
ing the BS equation is not new; 7 - we tread the 
subject in detail for the sake of completeness. 
The antiferromagnetic spin-dependent interaction 
J^2<i j> S-i-Sj = J\ + J2 consists of two terms: 

Ji = TE<i,j>[niAt + ™a%i - ^»,t«j4 - ^a«j,f] 



a position vector r^. The symbol X^<ij> means sum over an d J; 
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The interaction described by the J2 term does not 
allow us to use two-component Nambu fermion fields, 
and therefore, we introduce four-component Nambu 

fermion fields tp(y) = ^|(2/)V , J(2/)V , t(y)V , |(j/)J and 

tp(x) = ^|(x)^|(x)V't(^)V'|( :c )^ ; where x and y 
are composite variables and the field operators obey 
anticommutation relations. The "hat" symbol over any 
quantity O means that this quantity is a matrix. 

The interaction part of the extended Hubbard Hamil- 
tonian is quartic in the Grassmann fermion fields so the 
functional integrals cannot be evaluated exactly. How- 
ever, we can transform the quartic terms to a quadratic 
form by applying the Hubbard-Stratonovich transforma- 
tion for the electron operators: 8 



DAe 



[±A a {z)D^l-\z,z')Ae(z)+^(y)f^\y-, X \ Z )^(x)A a (z)\ 



±4>(y)ri HyMz)Hx)D^>(z,z')i,(y')r^>(y'-,x'\z')iP(x') 



(3) 



The last equation is used to define the 4x4 matri- 
ces D^j and fL 0) {a, (3 = 1,2,3,4). Their Fourier 
transforms, written in terms of the Pauli cr^, Dirac 
7 and alpha 7 7 matrices, are as follows: D^ ' — 



D 1 
D 2 



, tfl = {j°±a z )/2 and = (a x ±ia y )/2, 



(o) _ 



where a, = ( * ^ J , D, = ( J(k) - V(k)) a + 

(U - J(k) - V{k)) a x and D 2 = 2J(k)a x . For a 
square lattice and nearest-neighbor interactions V^(k) = 
4V(cos(k x ) + cos(ky)) and J(k) = J(cos(k x ) + cos(k y )). 



After applying transformation ([3]) the system under con- 
sideration consists of a four-component boson field A a {z) 

interacting with fermion fields ip(y) an d ip(x). The 



action of the model system is S 



_ c(e) 



5', 



S^ e ~ A ) where: sj\ e) = </> 



(A) 



tA„(z)D^-\z,z>)Ap(z') and = ^(y)f^(y,x \ 

z)ip{x)A a {z). Here, we have introduced composite vari- 
ables x,y,z — {r^it}, where is a lattice site vector, 
and variable u range from to j3 = \/kgP (T and ks 
are the temperature and the Boltzmann constant). We 
use the summation-integration convention: that repeated 
variables are summed up or integrated over. 

The idea of using the Hubbard-Stratonovich transfor- 
mation is to establish an one-to-one correspondence be- 
tween the t — U — V — J model and the model used 
to describe the Bose-Einstein condensate of excitons in 
semiconductors^. This allows us to follow the same steps 
as in Refs. j}], and to derive a set of 20 BS equations 
for the collective mode w(Q) at zero temperature. The 
Fourier transforms of V and J interactions are separa- 
ble, and therefore, we obtain a set of 20 coupled linear 
homogeneous equations for the dispersion of the collec- 
tive excitations. The existence of a non-trivial solution 
requires that the secular determinant cfei||x _1 — V\\ is 
equal to zero, where the bare mean-field-quasiparticle 

response function \ — ^ qT ^ ^ an d the interac- 
tion V = diag{U, U, -{U - 2J(Q)), U - 2F(Q), -{2V + 
J/2), -(2V+J/2), -{2V - J/2), -(2V - J/2)) are 
20 x 20 matrices. P and Q are 4x4 and 4 x 16 blocks, 
respectively, while R is 16 x 16 block (in what follows 
i,j = 1,2,3,4): 
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The quantities I a .b — F a ^(e(k, Q)) and J a ,b — F a ^(oj)^ the 1x4 matrices P a 



;, 6 (e(k,Q)) and J* b 



and the 4x4 matrices l %3 b = 



_^ 6 ( e (k,Q))and^ 6 
&(k, Q) = ^k,Q,m ki Q,7 k:Q or 7k,Cj): 



F % J b {uj) are defined as follows (the quantities a(k, Q) and 



F a ,b{x) = — 



xa(k,Q)b(k,Q) 



N ^ w 2 -e 2 (k,Q) 



F l 



M - 1 V- ™(k,Q)6(k,Q)A , k tJ 



TV 4^ 



xa(k,Q)6(k,Q) 
w 2 -e 2 (k,Q) 



Here e(k, Q) = E(k + Q) + E(k), and the form 
factors are: 7 k ,Q = Wk^k+Q + «k«k+Q, ^k,Q = 
Mk^k+Q - v k v k+ Q, 7k : Q = Ukt'k+Q - Mk+Q«k, and 
^k.Q = "k^k+Q + "k+Q^k where u 2 . = \ - v\ = 
[l + s(k)/E(k)]/2. 



We can compare our BS equations with the previ- 
ous approaches. If one takes into account only the spin 
channel, the secular determinant is a 1 x 1 determi- 
nant and the equation for the collective mode assumes 
the form 1 + (U - 2J(Q)7 77 = 0. This is the well- 
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FIG. 1: Lines in U, J parameter space which reproduce the 
resonance energy of 0.04 eV. The H&C line has been calcu- 
lated using the set of parameters given in Ref. [f| , while set 
1 and set 2 are obtained by using the two sets of parameters 
given in Ref. [ill ]. 



known equation for the poles of the spin susceptibility 
in the random phase approximation (RPA) x(Q; w ) = 
X (0) (Q;u;)/ [l + (U - 2J(Q)x (0) (Q;^)] , where the BCS 
susceptibility at zero temperature is x ™(Qi w ) = -£77^ 
Obviously, the RPA overestimates spin fluctuations be- 
cause the mixing between the spin channel and other 
channels is neglected. In the m— channel response the- 
ory we take into account the mixing between the spin 
channel and m — 1 other channels. The secular deter- 
minant is m x m. The coupling of the spin and two 
it channels leads to a three-channel response-function 
theory^ The four-channel theory^ includes the extended 
spin channel to the previous three channels. To see 
the difference between the previous theories and the 
present approach, we rewrite the secular determinant as 
A B 



detlx- 1 -V\ = det 



B C 



= det\C\det\A-BC~ l B q 



The 771— channel response-function theory takes into ac- 
count only the m x m matrix A, neglecting the mixing 
with the other 20 — m channels which is represented by 
the m x m matrix BC~ X B T . 

Strengths of interactions in BiTlYl samples. For 

Bi2212 compound, there are two possible sets of param- 
eters with all tight-binding basis functions involved (see 
Table 1 in Ref. [HI). Assuming A = 35 meV, we ob- 



tain Vd 1} = 115.2 meV with set 1, and VY' = 87.9 meV 
with set 2. Hao and Chubuko'v 6 - have used another set of 
parameters (we shall call it H&C) for Bi2212 compound 
with a doping concentration x = 0.12: t% = — 4i,i 2 = 
1.2t, t = 0.433 eV, ,u = -O.Mt, A = 35 meV and 
Y# — 0.6t. The parameters U, V and J should be ad- 
justed in such a way that the sharp collective mode of 
40 meV which appears at wave vector corresponds 
to the lowest collective mode calculated by the BS equa- 
tions. In Fig. 1 we present the results of our calculations 



of the lines in U, J parameter space which reproduce the 
resonance energy of 40 meV using all twenty channels but 
with set 1, set 2 and H&C parameters. Perhaps due to 
the limited size of single crystals currently available, no 
incommensurate peaks in Bi2212 samples have been re- 
ported so far. To obtain the exact values of the strengths 
we need at least one IC peak. 

Strengths of interactions in Y Ba^Cu-^O^.-j sam- 
ples. To obtain the strengths of the interactions we use 
the commensurate peak at ~ 40 meV and incommensu- 
rate peaks at ~ 24 meV and ~ 32 meV which have been 
reported in underdoped Y Ba-iCu-iO^.^. It is known that 
the YBaCuO is a two-layer material, but most of the 
peak structures associated with the neutron cross section 
can be captured by one layer band calculations^. The 
effects due to the two-layer structure can, in principle, 
be incorporated in our approach, but this will make the 
corresponding numerical calculations much more compli- 
cated. 

The commensurate and incommensurate peak struc- 
tures associated with the neutron cross section in 
YBaCuO have been studied within the electron-hole sce- 
nario using the single-band Hubbard t — U modeli^r— 
or the t — J modeler—. The techniques that have 
been used are based on (i) the Monte Carlo numeri- 
cal calculations^, (ii) the random phase approximation 
(RPA) for the magnetic susceptibility—^, (iii) the mean- 
field approximatio n 16 ' 17 and (iv) the RPA combined with 
the slave-boson mean field scheme^—. It is known 
that in the case when the Hubbard repulsion is large 
(U/t — > 00), the antiferromagnetic exchange J interac- 
tion is the consequence of Hubbard repulsion, because 
the t — J model is obtained after projecting out the dou- 
bly occupied states in the Hubbard t — U model, so that 
J = 2t 2 /U. Strictly speaking, by projecting out the dou- 
bly occupied states we remove the high-energy degrees of 
freedom and replace them with kinematical constraints 
assuming that the high energy scale (given by U that in 
cuprates corresponds to the energy cost to doubly occupy 
the same site) is irrelevant. Thus, if the constraint of no 
double occupancy is released, we arrive to the conclusion 
that the the magnetic susceptibility should be calculated 
using the t—U — J model rather than the t — U and t — J 
models (the corresponding arguments are presented in 

Refs. [ii-ii]). 

In our calculations the mean-field electron energy 
has a tight-binding form £k = — 2£ (cos k x + cos k y ) + 
At' cos k x cos k y — 2t"(cos2k x + COS2A;,,) — fi . Us- 
ing the established approximate parabolic relationship 
T c /T Ctmax = 1 - 82.6(p - 0.16) 2 , where Tc,max - 93K 
is the maximum transition temperature of the system, 
T c = 67K is the transition temperature for underdoped 
YBa-iCuzO^.T, we find that the hole doping is p — 0.10. 
At that level of doping the ARPES parameters are ob- 
tained in Refs. [H|: t = 0.25 eV, t' = 0.4/j, t" = 0.0444i 
and n = —0.27 eV. In the case of d-pairing the gap func- 
tion is Ak = Adk/2, where the gap maximum A should 
agree with ARPES experiments. In the case of under- 
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doped Y BaiCuzOs.i the gap maximum has to be be- 
tween the corresponding A = 66 meV in YBa2Cu^O^ ^ 
and A = 50 meV in Y BC12CU3O6.93H, so we set A = 60 
meV. The numerical solution of the gap equation pro- 
vides V[p — 265 meV. Next, we have solved numerically 
the BS equations to obtain the spectrum of the collec- 
tive modes w(Q) at the commensurate point Q^p, as 
well as at four incommensurate points 7r(l ± S),tt) and 
7r, 7r(l ± 6)). From the 40 meV solution of the BS equa- 
tions at Qa F we obtain a relation between U and J 
parameters, which is represented by the linear formula 
U = -3.985 J+1.01 (U and J are in eV). By means of the 
last relation we have solved the BS equations for the one 
of the four incommensurate 24 meV peaks (S — 0.22). 3 
The solution provides the following interaction strengths: 
J - 129 meV, V ~ 35.7 meV and U ~ 495 meV. 
To test the above values of the strengths we calculated 
the positions of the incommensurate peaks at 32 meV 
(5 = 0.19)4. The BS equations with the above strengths 
provide the deviation from of about 8 = 0.192, 

which is in excellent agreement with the experimentally 
obtained deviation (see FIG. 2 in Ref. Q). The strength 
of J is in a very good agreement with the strength of the 
superexchange interactions in the underdoped antiferro- 
magnetic insulator state of the cuprates which in YBaCO 
family has a magnitude of 0.1 — 0.12 cV, though some 
theoretical papers have predicted similar magnitudes. 

Strengths of interactions in La2-o.wSrQ,iQCuOi 
samples. To obtain the strengths of the interactions 
we use two IC peaks at ~ 7.1 meV (5 = 0.261) and at 
~ 10 meV (6 = 0.255) which have been reported in Refs. 
[28| . The tight-binding band parameters interpolated 
from values given in Refs. [Hj]: t = 0.25 eV, t' = 0.148*, 



t" = — 0.5i and ^ = — 0.821i. The maximum energy gap 
(estimated according to the prediction of the mean-field 
theory 3 ^) is about 10 meV, and therefore, Vy, = 95.8 
meV, which corresponds to maximum value of J about 
64 meV. The solutions of the BS equations provide the 
following interaction strengths: J ~ 44.1 meV, V ~ 14.8 
meV and U ~ 442.5 meV. To test these parameters we 
calculated the position of the IC peak at 18 meV. The 
calculated value of <5 = 0.19 is in very good agreement 
with the experimental value of S = 0.186*22, The above 
strengths provide the position of the commensurate peak 
to be at ~ 40 meV, which is in a very good agreement 
with the experimental value of 41 ±2.5 meV^ It is worth 
mentioning that in La^CuO^ J ~ 120 meV, and there- 
fore, our calculations support the idea that the antifer- 
romagnetic interaction decreases with doping. 31 
Conclusion Our approach to the incommensurate and 
commensurate structure of the magnetic susceptibility is 
based on the conventional idea of particle-hole excitations 
around the Fermi surface. There exists a second inter- 
pretation of the structure of the magnetic susceptibility 
in terms of the spin-charge stripe scenario (see Ref. [321 ] 
and the references therein) according to which the incom- 
mensurate peaks are natural descendants of the stripes, 
which are complex patterns formed by electrons confined 
to separate linear regions in the crystal. We do not wish 
to repeat the theoretical arguments that were advanced 
against the stripes model, but our unified description of 
the peaks based on the t — U — V — J model strongly 
supports the idea put forward by various groups that the 
commensurate resonance and the incommensurate peaks 
in cuprate compounds have a common origin. 
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